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Description 

[0001] L'invention concerne un dispositif et une methode pour g6nerer une fonction d'approximation, celle-ci etant 
fondee sur des premiers couples de valeurs associant une grandeur dependante k une grandeur independante, et 

5 pour determiner des seconds couples de valeurs desdites grandeurs a partir de ladite fonction d'approximation. 

[0002] Un dispositif et une methode de ce type sont connus du brevet US-A- 3 789 203 qui decrit un generateur de 
fonctions operant une approximation par interpolation iterative. Ce dispositif est prevu pour etre utilise dans des ap- 
plications de traitement de donnees necessitant un calcul de fonctions telles que sin(x), tg(x) par exemple. Ce dispositif 
ne requiert qu'une capacite de stockage minimale de la part d'un dispositif utilisateur. A partir de deux points appar- 

10 tenant k une fonction a interpoler, la methode tout d'abord interpole la fonction par une droite reliant les deux points, 
puis fait une approximation des ecarts entre la droite et la fonction par des approximations polynomials de degre 
croissant. Ensuite elle substitue aux points initiaux des points approximatifs pour reduire la longueur du segment reliant 
ies points k traiter et enfin r£itere les operations precedentes. 

[0003] Une telle methode necessite des ressources importantes en moyens de calcul et ne peut etre mise en oeuvre 
15 qu'avec des calculateurs performants. 

[0004] Or, il existe des applications ou une telle methode ne peut pas etre mise en oeuvre car elles ne disposent 
pas des ressources suffisantes. De plus, pour certaines applications on peut se satisfaire d'un calcui approche de la 
fonction pour des valeurs, en nombre limits, de la grandeur independante. 

[0005] II peut s'agir d'une fonction sigmoide appliquee a des potentiels neuronaux delivres par au moins un neurone 
20. dans un reseau de neurones. II peut s'agir d'une autre fonction non lineaire, par exemple une fonction racine, pour 
des calculs de distances entre des etats de neurones. Les applications peuvent aussi concerner d'autres dispositifs 
comme des generateurs de fonctions, des calculateurs ou autres. 

[0006] Pour calculer une telle fonction, sans passer par une fonction d'approximation, on peut utiliser differentes 
manieres. 

25 [0007] On peut eff ectuer le calcul mathematique exact pour chaque valeur de la grandeur independante a traiter, en 
programmant un calculateur selon les methodes connues. Une telle methode necessite d'effectuer k chaque fois les 
memes operations ce qui peut necessiter beaucoup de temps si le nombre de valeurs est eleve. 
[0008] On peut aussi pr6alablement stocker dans une memoire des tables precalcul6es. Dans ce cas, la lecture en 
memoire du resultat peut etre rapide. Mais pour couvrir avec un pas assez fin toutes les valeurs possibles de la grandeur 

30 independante, il faut alors disposer de tables de grandes capacites. Ces methodes de calcul presentent done des 
inconv6nients. 

[0009] D'autre part, on peut etre conduit k identifier deux grandeurs qui sont dans la dependance Tune de I'autre par 
des couples de valeurs associant une grandeur dependante k une grandeur independante. Ainsi, dans le suivi d'un 
processus industriel, on peut etre conduit k mesurer par exemple un rendement R d'une operation en fonction de la 

35 temperature T a laquelle a ete realisee ladite opeation R = f(T) . Un suivi du processus peut consister k consigner des 
lots de couples de mesures sur un graphique. Celui-ci peut etre utilise pour caracteriser le processus ou pour en 
deduire de nouveaux parametres de reglage de ladite operation. Ceci fait par exemple I'objet de I'article de H. ISHI- 
BUCHI et H. TANAKA, "Regression analysis with interval model by neural networks" in "IEEE International Joint Con- 
ference on Neural Networks", vol. 2, 18-21 Nov. 1991, SINGAPORE. II est souhaitable que ces nouveaux parametres 

40 soient representatifs du fondement de ladite operation et que les fluctuations de mesures inherentes k ce genre de 
processus soient 6cart6es. II est done souhaitable de determiner une approximation de la fonction f(.). 
[001 0] Ainsi dans un cas il peut s'agir de mesures erratiques ou entachees d'erreurs que Ton desire representer par 
une fonction d'approximation. 

[0011] Dans un autre cas, on connait des valeurs precises mais ('utilisation ci en faire ne necessite pas une grande 

45 precision et une fonction d'approximation suffit. 

[0012] Un des buts de l'invention est de generer une fonction d'approximation avec des moyens materiels r6duits 
permettant de calculer rapidement un nombre reduit de valeurs de la grandeur dependante utiles k ('application sans 
avoir pour cela k determiner d'autres valeurs de la fonction d'approximation. Un but compl6mentaire est de d6livrer 
des valeurs qui peuvent etre approchees dans la limite d'une erreur maxima le controlee. 

50 [0013] Ce but est atteint avec un dispositif caracterise" en ce qu'il comprend : 

des premiers moyens : 

pour determiner iterativemant au moins une fonction lineaire courante de regression en rendant egales, en 
55 valeur absolue, des premieres erreurs de signes alternes mesur£es entre des premieres valeurs de la grandeur 

dependante pour trois couples d'une suite desdits premiers couples, et respectivement des secondes valeurs 
de la grandeur dependante determinees, d'apres ladite fonction lineaire courante, pour les memes valeurs de 
la grandeur independante, 
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pour selectionner celle des fonctions lineaires courantes qui deiivre ('approximation de tous les couples de 
ladite suite avec des erreurs minimales, 

el pour coder, a I'aide de codes specifiques, la fonction lineaire de regression selectionn6e, 

s - et des seconds moyens pour determiner lesdits seconds couples a I'aide desdits codes specif iques. 

[0014] Ainsi avantageusement on determine une fonction lineaire de regression approchant au mieux les differents 
couples de valeurs connus. Les resultats approch6s ainsi deiivr6s forment un compromis satisfaisant pour de nom- 
breuses utilisations du dispositif g6n6rateur de fonction. 
w [0015] Une fonction lineaire de regression est une fonction simplificatrice qui represente un phenomene complexe 
en reduisant les parametres significatifs. En reprSsentant la suite de couples de valeurs par des points dans un espace 
a deux dimensions, la fonction lineaire de regression devient une droite de regression. 

[0016] Ainsi apres avoir defini la droite de regression par des codes, on peut calculer une valeur approch6e de la 
grandeur dependante en tout point de la droite de regression avec des moyens reduits pour des valeurs quelconques 
75 de la grandeur independante. 

[0017] L'invention concerne egalement une methode pour genSrer une fonction d'approximation , la methode 
comprenant : 

une premiere phase : 

20 

pour determiner iterativement au moins une fonction lineaire courante de regression en rendant egales, en 
valeur absolue, des premieres erreurs, de signes alternes, mesurees entre des premieres valeurs de la gran- 
deur d6pendante pour trois couples d'une suite desdits premiers couples, et respect ivement des secondes 
valeurs de la grandeur d6pendante determinees, d'apres ladite fonction lineaire, pour les m£mes valeurs de 
25 la grandeur independante, 

pour selectionner celle des fonctions lineaires courantes qui deiivre I'approximation de tous les couples de 
ladite suite avec des erreurs minimales, 

et pour coder la fonction lineaire de regression s6lectionn6e a I'aide de codes specifiques, 

30 - et une seconde phase pour determiner lesdits seconds couples a I'aide desdits codes specifiques. 

[001 8] Les moyens mis en oeuvre par l'invention peuvent etre formes par un calculateur programme ou par un circuit 
d6di6. Ms peuvent aussi mettre en oeuvre des neurones. 

[0019] Un dispositif mettant en oeuvre des neurones selon l'invention peut etre utilise par un reseau de neurones, 
35 dont il peut notamment en constituer un sous-ensemble. En effet, pour fonctionner, le reseau de neurones doit disposer 
de moyens pour appliquer une fonction non lineaire d'activation aux potentiels de neurones qu'il deiivre. Selon l'inven- 
tion, le dispositif muni de neurones peut calculer une approximation de cette fonction non-lin6aire d'activation. II peut 
6galement calculer des distances entre des 6tats de neurones en calculant une approximation d'une fonction racine 
carree destinee a §tre exploitee dans le reseau de neurones. 
40 [0020] Lorsque la taille de la suite de couples de valeurs fournis initialement est elevee, on peut diviser la suite de 
couples en plusieurs sous-ensembles pour determiner plusieurs droites de regression et ameiiorer la precision de 
I'approximation. La fonction d'approximation de la suite de couples est alors formee par une fonction lineaire par mor- 
ceaux pour laquelle une exigeance de continuite entre les morceaux peut etre ou non imposee. 
[0021] Certains couples de I'ensemble de couples de valeurs peuvent avoir une influence particuliere que Ton peut 
45 concretiser en donnant un coefficient de ponde ration specifique a chaque couple. Dans ce cas, I'erreur affectee a 
chaque couple tient compte de ce coefficient de ponderation specifique. 

[0022] Ces differents aspects de l'invention et d'autres encore seront apparents et eiucid6s a partir des modes de 
realisation d6crits ci-apres. L'invention sera r6alis6e selon les revendications independantes 1 et 6 ci-jointes. 
[0023] L'invention sera mieux comprise a I'aide des figures suivantes donnees a titre d'exemples non limitatifs qui 
50 represented : 

Figure 1 : un graphique montrant une representation a deux dimensions d'un ensemble de points avec une droite 
de regression D. 

Figure 2 : un graphique montrant un ensemble de points et des droites servant a la determination d'une enveloppe. 
55 Figure 3 : un organigramme d'une premiere variante de mise en oeuvre de la methode a partir de triplets de points. 

Figure 4 : un organigramme d'une seconde variante de mise en oeuvre de la methode a partir de triplets de points. 
Figure 5 : une partie d'organigramme d'une troisieme variante de mise en oeuvre de la methode a partir de couples 
de points. 
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Figure 6 : une partie d' organigramme indiquant une preselection de points appartenant a une enveloppe inferieure 
ou a une enveloppe superieure de I'ensemble des points. 

Figure 7 : un organigramme pour la determination des enveloppes inferieure et superieure. 
Figure 8 : un schema d'un dispositif selon I'invention. 
s Figure 9 : un schema d'un dispositif de codage pour le calcul et le codage d'une droite de regression dans.le cas 

general. 

Figure 10 : un schema d'un dispositif de codage pour le calcul et le codage d'une droite de regression pour la 
premiere variante. 

Figure 1 1 : un schema d'un dispositif de transcodage pour le calcul des valeurs de la grandeur d6pendante a partir 
10 d'un ensemble de droites de regression cod6es. 

Figure 12 : un schema d'une realisation neuronale des moyens de comparaison de la figure 10. 

Figure 13 : une representation montrant une fonction d'approximation formee de plusieurs droites. 

Figure 14 : deux representations montrant une determination des limites de validite de deux droites consecutives. 

Figure 15 : une representation concernant le raccordement de droites consecutives. 

15 

[0024] L' invention concerne I'approximation d'une fonction connue uniquement a travers un certain nombre de points 
par exemple P v .. P 6 (figure 1) dans une representation a deux dimensions. Chacun de ces points est defini par un 
couple de valeurs (x, y) reliant la gfandeur independante x a la grandeur dependante y. Par la suite, il sera question 
d'une paire (respectivement d'un triplet) de points, ce qui met en oeuvre deux couples (respectivement trois couples) 
20 de valeurs. On ordonne les points d'apres un ordre croissant des valeurs d'abscisses Xj, ce qui definit un indice i 
croissant avec lesdites valeurs. Une convention inverse peut etre faite en modifiant en consequence Texpose ci-apres. 
[0025] Selon I'invention on effectue une approximation de I'ensemble des couples (X^ Y-, ), (X 2 , Y 2 )... par une droite 
de regression D ayant pour equation : 

25 

(1) D:y=p.x+q. 

ou x et y sont des variables courantes. 

[0026] Pour cela on considere dans une premiere variante trois couples de valeurs, par exemple (Xq, Y 3 ), (X^ Y 4 ), 
30 Y 5 ), et on determine une droite de regression D par equilibrage des erreurs absolues. Une erreur est mesuree 

par la difference apparaissant, pour une abscisse x donnee, entre la valeur y du point et I'abscisse y mesuree sur la 
droite de regression. Equilibrer les erreurs sur trois points, consiste a avoir trois erreurs egales en valeur absolue avec 
un signe d'erreur oppose aux deux autres pour le point ayant une abscisse x comprise entre les abscisses x des deux 
autres points. Puis on examine si pour les points restants de I'ensemble, Perreur qui les separe de la droite reste 
35 inferieure ou egale en valeur absolue a I'erreur pr6alablement determin6e pour les trois points seiectionnes. Ceci 
decoule du fait que Ton s'interesse a une erreur pire cas relative a I'ensemble de tous les points a prendre en consi- 
deration c'est-a-dire la plus grande erreur, en valeur absolue, qui existe entre un des points et la droite de regression. 
Si toutes les erreurs sont effectivement inf6rieures ou egales, la droite est s6lectionn6e pour repr6senter les points 
sinon on recommence les operations avec trois nouveaux couples de valeurs pour determiner une autre droite de 
<o regression. 

[0027] II peut exister plusieurs droites de regression representant tous les points de I'ensemble. Selon (a methode 
on determine la droite de regression optimale qui minimise I'erreur pire cas definie prealablement. 
[0028] La figure 1 represente un exemple forme de six points P t a P 6 disposes selon une representation a deux 
dimensions. A des fins duplication considerons ie resultat final. On observe que la droite de regression D de la figure 

45 1 est situee de sorte que les erreurs sont egales en valeur absolue pour les points P 3 , P 4 et P 5 . Pour les points P v 
P2. P6 les erreurs sont inferieures aux precedentes, en valeur absolue. Dans un ensemble de points P, a P 6 , la methode 
va ainsi consister a rechercher les trois points particuliers, ici P 3 , P 4 , P 5 , qui permettent de determiner la droite de 
regression optimale qui minimise I'erreur pire cas puis a coder cette droite. Dans le cas de la figure 1 representant un 
resultat final, si on trace deux droites D-, et D 2 , paralleles a la droite de regression D, qui passent respectivement par 

50 les points P 3 , P 5 , d'une part et P 4 d'autre part, on constate que tous les points de I'ensemble sont a I'interieur d'un 
bandeau limite par les droites D 1 et D 2 ou sur ces droites. 

[0029] La phase de determination de la droite de regression peut donner lieu a plusieurs mises en oeuvre dont seules 
les plus avantageuses seront decrites ci-apres. 

55 Premiere variante de la premiere phase de la methode . 

[0030] La figure 3 represente la suite des etapes a mettre en oeuvre pour determiner la droite de regression optimale. 
[0031] Parmi I'ensemble des points P v .. P N (bloc 100), on seiectionne (bloc 102) trois points quelconques P v P;, P k 
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avec i < j < k. Ces trois points servent a determiner la droite de regression qui minimise I'erreur sur y pour ces trois 
points. Cette determination est effectuee de maniere analytique preferentiellement par des moyens programmes. On 
determine une droite de regression D pour que trois erreurs (bloc 104) : 

E PD (P l ,D),Ep D (P j ,D) ( Ep D (Pk > D) 
entre la droite D et chacun des points verifient : 

E PD (P i( D) = -E pD (P., D) = E pD (P k .D) 

avec 

E PD (P jt D) = Y j -(p.X i+ q) 
et des relations analogues pour les autres erreurs. 

[0032] On determine la droite de regression D a J'aide des coefficients p et q de Tequation (1) tels que : 

_ V Yj 

P ~ Y - Y 

et 

X y . Y f - X r Yj + X k Yj - Xj. Y k + X^ Yj - X- r Y k 
q - 2(X k -X.) 

[0033] L'erreur se rapportant a un triplet (Pj, P|, P k ) s'ecrit alors : 

E T (P,Pj,P k ) = \E pD ( Pi ,Dl 

[0034] Lorsque I'erreur E T a 6te ainsi calcul6e, on examine si les autres points de I' ensemble engendrent des erreurs 
inf 6rieures ou egales, en valeur absolue, h celles des points P j( Pj, P k . Pour cela on s^lectionne un point P m (bloc 1 06) 
additionnel et on calcule la valeur absolue de I'erreur E Pm (bloc 108) apparaissant entre la valeur de la variable y au 
point P m et la droite D. 

[0035] Lorsque cette erreur E Pm est inferieure ou 6gale en valeur absolue a E T (bloc 110) (repere Y), le point P m 
additionnel est accepte et la methode se poursuit (bloc 112) avec un point additionnel suivant (bloc 106). Si tous les 
points additionnels satisfont au critere I I < E T , la droite D est acceptee et ses coefficients sont utilises pour coder 
la droite de regression optimale D opt (bloc 114). 

[0036] Lorsque cette erreur E Pm est superieure a E T (bloc 110) (repere N), le triplet de points Pj, Pj, P k selectionne 
n'est pas accepts et un autre triplet de points (bloc 116) est choisi (lien 101) dans I'ensemble de points (bloc 102). La 
methode s'acheve & I'obtention de la droite D (bloc 114) satistaisant ce critere meme si tous les triplets de points n'ont 
pas ete examines. 

[0037] 1 1 peut apparaitre des situations ou & Tissue de I'etape 1 1 6 tous les points possibles ont ete examines et aucun 
des triplets n'a fourni de solution (bloc 118). Dans ce cas, il est possible de reprendre le deroulement de la premiere 
variante en remplacant le test du bloc 110 par le test suivant : 

ou a est un coefficient legerement superieur a 1 . Dans cette meme situation, il est egalement possible de faire appel 
a la deuxieme variante. 

[0038] La premiere phase peut comprendre les etapes suivantes : 
A - selection de trois couples de valeurs parmi ladite suite, 

B - calcul de la lonction lineaire courante de regression D et determination d'une erreur de triplet E^IEp^l s'y 
rapportant, 
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C - selection d'un couple additionnel, 

D - calcul d'une erreur additionnelle E Pm entre le couple additionnel et ladite fonction, 

E - et lorsque I E Pm I < £^1 1 0) pour le couple additionnel, la methode reprend a I'etape C avec un couple additionnel 
suivant, 

F - et lorsque I E Pm I > E T pour au moins un couple additionnel, 

la methode reprend a I'etape A avec une selection d'un autre groupe de trois couples parmi ladite suite de couples, 

G - et lorsque I E Pm I < E 7 pour tous les couples additionnels, la fonction lineaire courante de regression est codee 
et stockee en tant que fonction lineaire ^approximation. 

[0039] On peut choisir de scruter P.ensemble des triplets en prenant un ordre croissant ou un ordre decroissant ou 
un ordre aleatoire pour effectuer cet examen. Le triplet qui sera retenu pour determiner la droite de regression pourra 
de ce fait etre detecte a un instant quelconque du deroulement de cette scrutation. II s'ensuit que la rapidite d'obtention 
de la droite de regression depend de I'instant au cours duquel le triplet est detecte. Sa mise en oeuvre presente un 
degre de complexite allant de N a N 4 ou N est le nombre de points initiaux. Sa complexity est done reduite pour un 
petit nombre de points. Cette variante permet d'obtenir une realisation materielle avec une forte parallelisation. Elle 
est tres peu sensible a une troncature des valeurs et fournit un resultat exact. 

Deuxieme variante de (a premiere phase de la methode. 

[0040] Dans cette deuxieme variante (figure 4), on selectionne successivement des triplets de points, on calcule 
chaque fois une droite de regression, et, par recurrence, on selectionne celle qui delivre Perreur E PD la plus grande 
e'est-a-dire correspondant a Perreur pire cas pour Pensemble de points consideres. 
[0041] Selon la seconde variante, la premiere phase comprend les etapes suivantes : 

A - selection de trois couples de valeurs parmi ladite suite, 

B - calcul d'une fonction lineaire courante de regression D et determination d'une erreur de triplet E-^\E pD \ s*y 
rapportant, 

C - comparaison de Perreur E T avec une erreur optimale E op initialisee a une valeur strictement negative, 

D - et lorsque E T > E op , mise a jour de Perreur optimale E op en remplacant E op par E T et mise a jour des codes 

d'une fonction lineaire optimale de regression D op en remplacant ceux-ci par les codes de la fonction lineaire 

courante D, 

E - puis retour a I'etape A pour selectionner trois autres couples, 

F - et lorsque tous les triplets de couples de valeurs de la suite ont ete testes, les demiers codes de la fonction 
lineaire optimale D op constituent les codes de la fonction lineaire d'approximation. 

[0042] A chaque examen d'un triplet, on compare Perreur E T dudit triplet avec Perreur optimale precedemment me* 
morisee et on met a jour Perreur optimale E op avec la plus grande valeur de Perreur E T determinee pour chaque triplet. 
On met a jour egalement les parametres de ia droite optimale. Avant la mise en oeuvre de I'etape A, il taut initialiser 
la valeur E op a une valeur faible negative, par exemple -1 . 

[0043] Dans ce cas, la rapidite d'obtention de la droite de regression ne depend pas du mode de scrutation des 
triplets. Sa complexite est de degre N 3 . Une mise en oeuvre materielle peut beneficier de la grande regularity de 
Palgorithme mis en oeuvre. II est peu sensible aux troncatures des donnees et fournit une solution exacte. 

Troisieme variante de la premiere phase de la methode. 

[0044] Dans cette troisieme variante (figure 5), on selectionne d'abord une pa ire de points a laquelle on ajoute un 
point supplemental, situe entre ces deux points, afin de former un triplet de points. Pour cela, on modifie les etapes 
A, B et C de la premiere variante, les autres etapes restant les memes. Les etapes modifiees sont telles que : 

A1 - modifie I'etape A en operant une selection de deux couples de valeurs appartenant a ladite suite, tel qu'il 
existe au moins un couple additionnel intermediate ayant une grandeur independante (X 1 - Xg) comprise 
entre les grandeurs independantes dudit couple pour constituer au moins un triplet de couples, 

A2 - modifie I'etape A premierement en determinant une fonction lineaire annexe qui contient les deux couples 
selectionnes et deuxiemement en determinant des secondes erreurs entre les grandeurs dependantes des 
couples intermediaires possibles et ladite fonction lineaire annexe : 
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et lorsque ces secondes erreurs sont toutes de meme signe, selection du couple intermediate fournissant 
la plus grande seconde erreur, en valeur absolue, pour lormer un triplet de couples de valeurs forme du 
couple intermediaire et des deux couples selectionn^s, 

et lorsque ces secondes erreurs sont de signes diff6rents, reprise de la methode a I'etape A1 , 

5 

B1 - I'etape B est effectuee avec ledit triplet selectionne, 

C1 - modifie I'etape C en seiectionnant un couple additionnel dont la grandeur independante n'est pas comprise 
entre la grandeur independante des deux couples sSlectionnes. 

10 [0045] Lorsque I'erreur I E Pm I est superieure a I'erreur E T (bloc 110), la methode reprend a I'etape A1 (lien 101) 
avec une nouvelle selection de paire de points (bloc 102a). 

[0046] On observe que le deroulement de cette troisieme variante depend de la scrutation des valeurs et done des 
valeurs elles-memes. La complexity de la mise en oeuvre de cette variante varie entre N et N 3 ce qui lui confere un 
certain avantage par rapport aux variantes precedentes. Une forte parallelisation des moyens de mise en oeuvre peut 
75 etre operee mais Timplementation des moyens peut presenter un manque de regularity ce qui peut constituer un 
handicap pour une realisation integr6e. Cette variante est peu sensible aux erreurs d'arrondis des valeurs et fournit 
une solution exacte. 

Quatrieme variante de la premiere phase de la methode. 

20 

[0047] Elle concerne la determination de la droite de regression a partir des enveloppes. 

[0048] II est possible de reduire le nombre de triplets a examiner en determinant des enveloppes respectivement 
superieure et inferieure entourant les points extremes dans la representation bidimensionnelle de Pensemble de points. 
Une enveloppe superieure ou une enveloppe inferieure est definie telle qu'en joignant par une droite deux points 

25 adjacents quelconques de I'enveloppe, tous les autres points soient situes d'un meme cote respectivement de I'enve- 
loppe superieure ou de I'enveloppe inferieure. On determine ainsi tous les points appartenant a ces dites enveloppes. 
[0049] La determination de la droite de regression va consister a considerer les paires de points adjacents d'une 
des enveloppes auxquels on associe un point intermediate n'appartenant pas a ladite enveloppe pour constituer un 
triplet et op6rer comme cela vient d'etre decrit dans le cas des paires de points de la troisieme variante. Si une solution 

30 optimale n'a pas ete trouvee, on considere les paires de points de Pautre enveloppe. 

[0050] Pour mettre en oeuvre une enveloppe, on selectionne une paire de points adjacents appartenant a I'enve- 
loppe. On determine alors s'il existe un point dispose de telle facon que son abscisse soit intermediaire entre les 
abscisses des points seiectionn6s. Lorsque ce point n'existe pas on passe a une autre paire de points de la meme 
enveloppe. Pour certaines paires, lorsqu'il apparait qu'il existe un ou plusieurs de ces points intermediates, on choisit 

35 le point intermediaire le plus eioigne de la droite contenant les deux points de la paire pour former un triplet et pour 
determiner une droite de regression. Pour determiner si cette droite de regression peut etre seiectionnee comme droite 
de regression optimale, la methode met en oeuvre les mdmes operations que celles decrites pr6alablement dans le 
cas de la troisieme variante. 

[0051] Pour cela on modifie la troisieme variante telle que (figure 6), pr6alablement a I'etape 102a (figure 5), la 
*o premiere phase de la methode comprend une etape (bloc 100a) de determination d'une enveloppe inferieure et/ou 
d'une enveloppe superieure r6unissant les points les plus extremes de ('ensemble de points, la selection des paires 
de points a I'etape 102a etant faite parmi les points adjacents- appartenant a I'une ou I'autre enveloppe. La selection 
de ladite paire de points est effectuee lorsqu'il existe au moins un point intermediaire ayant une abscisse situ6e entre 
les abscisses des points de la paire de points. S'il existe plusieurs points intermediaires, on forme le triplet avec le 
<5 point intermediaire le plus eioigne de la droite passant par les deux points qui torment la paire de points. Si une solution 
n'est pas trouvee avec la premiere enveloppe, on poursuit le traitement avec la seconde enveloppe. 
[0052] La determination des enveloppes inferieure et superieure est effectuee selon I'organigramme de la figure 7. 
L'indice des points allant croissant avec la variable x, le premier point P c fait ainsi partie des deux enveloppes. Con- 
sid6rons d'abord I'enveloppe inferieure, les points appartenant a I'enveloppe inferieure sont rep6res par la lettre Q. Un 
50 point courant Q v est repere par l'indice v. 

[0053] A I'etape 400 les deux premiers points de I'enveloppe sont : Q Q = P 0 et Q 1 = P 1 . Le dernier point courant 
Q v est repere : v = 1 pour P v Le nombre kde points traites est comptabilise. 

[0054] Etape 402 : un premier test est effectue pour determiner si le dernier point P N . n a 6te traite pour detecter la 
fin de la determination de I'enveloppe. 
55 [0055] Etape 404 : dans le cas contraire on teste si v > 1 . Lorsque v < 1 on incremente v (v = v + 1 ) et on prend le 
point P k comme point Q v (etape 407). On incremente l'indice k pour traiter le point P suivant (etape 409) et reprise de 
la methode a I'etape 402. 

[0056] Lorsque v ^ 1 , on calcule une droite passant par les points . 1 et Q v (etape 406) et on determine le signe 
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e de I'erreur sur la variable dependante entre le point P k courant et cette droite (Q v 0 , Q v ). Ceci a pour but de determiner 
si le point courant est au-dessus ou en dessous de la droite (Qy . 1f Q v ). 

[0057] Lorsque le signe e < o, il faut supprimer le dernier point C^, decr6menter v tel que v = v - 1 (6tape 410) et 
reprendre la m6thode a l'6tape 404. On peut de cette maniere etre conduit a supprimer certains points deja accepter 
lorsqu'un point suivant necessite de !es supprimer. 

[0058] Lorsque le signe e est strictement positif, la methode rep rend a I'etape 407 avec un point suivant. 

[0059] Cet organigramme s'applique a la determination des enveloppes inferieure et superieure en inversant le signe 

de I'erreur a considerer. 

[0060] Pour faire comprendre les mecanismes ainsi mis en oeuvre considSrons, a tit re d'exemple, le cas simple 
forme par les points P v P 2 , P 3 , P 4 de la figure 2 et determinons I'enveloppe inferieure. Le point P, est le premier point 
de I'enveloppe d'ou Q 0 = P 1 . On prend P 2 comme deuxieme point . On calcule la droite L Q reliant P t , P 2 . Le segment 
P.,, P 2 est accepte comme premier segment de I'enveloppe basse. On examine si P 3 est au-dessus de la droite L 0 , ce 
qui est le cas sur la figure 2. Le segment P 2 P3 est accepte comme deuxieme segment de I'enveloppe basse. On 
calcule ensuite la droite L A reliant P 2 et P 3 et on examine si P 4 est au-dessus de la droite L v Ce n'est pas le cas sur 
la figure 2. Ceci entrame que le point P 3 et la droite sont d'abord ecart6s. De meme, P 4 est en dessous de la droite 
L Q done le point P 2 et la droite L c sont ensuite aussi ecartes comme ne pouvant pas appartenir a I'enveloppe inferieure 
et les operations se poursuivent avec la determination de la droite reliant P v P 4 , puis Pexamen des points suivants 
(non represents). 

[0061] La complexite de mise en oeuvre de cette variante basee sur des enveloppes est de degre N 2 et depend de 
I'ordre de scrutation des donnees. Cette complexite est moindre que celle des variantes precedentes et de ce fait 
fournit un resultat rapidement. La regularity de I'implementation est moyenne mais cette variante est peu sensible aux 
erreurs d'arrondi des valeurs et delivre une solution exacte. 

[0062] Pour certaines applications, il peut etre souhahable d'accroTtre la precision de la determination de la.fonction 
d'approximation dans certains domaines de la grandeur independante x et d'affecter des coefficients de pond6ration 
Wj aux points P r C'est par exemple le cas lorsque la fonction d'approximation est faiblement variable avec la grandeur 
independante x. A certains points peuvent alors etre affectes des coefficients de ponderation. Ceux-ci peuvent etre 
communs a plusieurs points ou etre individuels pour chaque point. Ces coefficients de ponderation Wj sont, pafr la 
suite, consideres strictement positifs. 

[0063] Dans les cas ou des coefficients de ponderation existent, on definit une erreur E PD entre un point Pj et la 
droite de regression D tel que : 

E ro (^D)=lV,[y r (p.X. + g)] 

ou E PD est une vaieur signee. La determination de la droite de regression D pour trois points Pj, Pj, P k est alors modifiee 
en ce que les quantity p et q de r equation (1 ) deviennent : 

(2) p = NUMP/DET et q = NUMQ/DET 
ou les quantites NUMP, NUMQ et DET sont definies par : 

DET = Wj.W.^X.-X.) + Wj.W^fX^Xj) + W-.W^fX^Xj) 

NUMP = Wj.Wj.fYj-Yj) + W r W k .(Y k -Y.) + W j .W k .(Y k -Y j ) 

NUMQ = W j .W j .(X i .Y r X i .Y j ) + W j .W k .(X k .Y r X j .Y k ) + 
W,W k .(X k .Y.-X,Y k ) 

[0064] Par ailleurs, I'erreur E T associee a ce triplet peut etre exprimee et calculee par : 
avec 
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E 1 iP i ,P f P k ) = \E PD lP i ,q 



E PD (Pi,D) = -E PD (Pj,D) = EpoiPpD) 



on 



(J) fc r= DET 

[0065] La premiere et la seconde variante de la methode de la premiere phase d6crites prScedemment peuvent §tre 
mises en oeuvre en faisant intervenir les coefficients de ponderation ci-dessus. Cette mise en oeuvre peut §tre effective 
is en programmant un calculateur. 

[0066] Dans le cas ou Ton seiectionne une paire de points en faisant intervenir des coefficients de ponderation W, 
la methode presente les amenagements suivants. 

[0067] Pour les points P jf P,, P k qui torment le triplet, on definit un facteur F, qui lie les coordonnees respectives et 
les coefficients de ponderation respectifs des points P { , P,, P k . Au point central 1 , on associe le facteur F, tel que : 

20 




25 

[0068] Ce facteur F, intervient dans la determination du point intermediate P, qui sera a seiectionner (i < I < k) pour 
constituer le triplet en tenant compte des poids affectant les points. On choisit un point et on calcule la droite de 
regression associee au triplet Pj, Pj, P k et I'erreur Ej| associee au triplet. 

[0069] Cette variante de la methode consiste a modifier uniquement l'6tape A2 de la troisieme variante (bloc 1 02b, 
30 figure 5). Cette etape determine I'existence et la valeur d'un point intermediate servant a former un triplet. On cherche 
d'abord a former une droite de regression situ6e en dessous de Pj et de P k . Pour chaque point intermediaire P () on 
determine si F, = 1 et E PD (Pj, D,) < 0. Si au moins un point verifie cette condition, ii n'existe pas de droite de regression 
situee en dessous de Pj et P k> sinon on determine une grandeur G max qui est la quantity maximale parmi : 

35 . d'une part les quantites E PD (Pj, D p ) pour tous les points intermediates, 

d'autre part les quantites (a) suivantes, uniquement pour les points intermediates pour lesquels F, < 1 avec 

40 h ' 1 

[0070] S'il existe au moins un point intermediaire tel que F ( > 1, on definit egalement une grandeur G min qui est la 
valeur minimale prise par (a) uniquement pour les points intermediaires tels que F| > 1. On examine alors s'il existe 
au moins un point intermediaire tel que : 
45 Ep D (Pj,D,)^0 

©max £ « E TI 

et G min . a > Ej, (a, coefficient £ 1). S'il existe un tel point, il est choisi comme point intermediaire pour former le 
triplet. 

[0071] S'il n'existe pas de point intermediaire tel que F-, > t , on teste s'il existe au moins un point intermediaire tel que : 
so Ep D (Pj,D,)£0 

G max - a E n ( a : coefficient £ 1 ). 
[0072] Si un tel point existe, il est choisi comme point intermediaire pour former le triplet. 

[0073] Si aucun triplet n'a et6 constitue, on cherche a former une droite de regression situee au-dessus des points 
Pj et P k . La m§me methode est reprise en inversant le signe des erreurs E PD . 
55 [0074] Si aucun point P, ne peut etre seiectionne, on recommence avec une autre paire P h P k . 

[0075] Lorsque I'ensemble de points a trailer est trop important pour etre represente par une seule droite de regres- 
sion, la methode code alors plusieurs droites de regression chacune etant determinee selon la methode d6crite pre- 
cedemment. 
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[0076] La figure 1 3 represente un exemple dans lequel la fonction d'approximation est constitute de plusieurs droites 
de regression. 

[0077] Dans une premiere situation, de par la connaissance des grandeurs a traiter, on peut vouloir imposer des 
limites a chaque droite de regression suivant la grandeur independante x. Ainsi on peut vouloir disposer d'une droite 
5 de regression D a entre les valeurs [X^XJ de la grandeur, borne x a incluse, borne x b exclue. De meme avec D b pour 
[X&, XJ et D c pour [X c , XJ . Dans ce cas le probleme revient a determiner une droite dans un domaine limits et a 
appliquer a chaque fois la methode deja decrite. 

[0078] Mais il est possible de faire que pour chaque droite, la methode determine des limites suivant la grandeur 
independante x sans qu'elles soient fixees au depart. 

w [0079] Le principe de la determination d'une limite optimale entre deux droites de regression adjacentes est repre- 
sente sur la figure 14-A. Soient deux droites de regression D1 et D2 non optimales. La droite D1 est determinee a 
partir de N1 points et la droite D2 est determinee a partir des N2 points restants avec N1 + N2 = N , ou N est le nombre 
total de points. Les droites D1 et D2 torment une approximation avec respectivement des erreurs maximales E1 et 
E2. Supposons que E2 soit superieure a E1 et que Ton cherche a reduire E2, i) faut alors transferer vers D1 un point 

is appartenant a D2. Lorsque le nombre de points diminue, I'erreur resultante est soit inchangee soit reduite. Lorsque le 
nombre de points augmente, I'erreur resultante est soit inchangee soit augmentee. En fonction du nombre n de points 
appartenant aux droites D1 et D2, les variations des erreurs E1 et E2 peuvent presenter les allures representees sur 
la figure 14-B pour I'exemple choisi. II s'ensuit que I'approximation obtenue globalement par D1 et D2 est optimale 
lorsque ces erreurs sont voisines, c'est-a-dire situees dans la zone deiimitee par un trait discontinu sur la figure 14-B. 

20 Ces erreurs ne sont pas necessairement egales compte tenu du caractere discontinu des valeurs de la grandeur 
independante. 

[0080] Pour determiner la valeur d'abscisse Xj im entre les deux droites : 

on determine la droite D1 sur un certain nombre de points et on calcule I'erreur E1 maximal©, 
2$ - on determine la droite D2 sur les points restants et on calcule I'erreur E2 maximale, 

on compare E 1 et E2 et on transffcre un point de la droite qui presente la plus forte erreur vers la droite qui presente 
la plus faible erreur, 

on determine la valeur limite X| jm lorsqu'il se produit une inversion dans ie rapport entre les deux erreurs. 

30 [0081] Cette methode s'applique par recurrence a plusieurs droites deiimitees entre elles par plusieurs limites d'abs- 
cisses. 

[0082] Les droites sont d6termin6es a partir d'une suite discrete et limitee de grandeurs de mesure. Neanmoins, 
pour ['exploitation des droites de regression, il est n6cessaire de definir leur domaine d'existence qui s'6tend sur un 
continum de valeurs situees entre deux limites d'abscisses. Or la determination des droites fournit une suite de droites 

35 qui ne se raccordent pas necessairement par leurs extremites. II peut apparaitre utile pour certaines applications 
d'eviter qu'un ecart apparaisse sur la valeur de la grandeur dependante y pour des valeurs voisines (Xp m - e) et (X, jm 
+ e) de la grandeur independante ou e est une tres faible valeur. II est possible de choisir que Pabscisse limite X| jm 
appartienne exclusivement a Tune ou a I'autre droite. II est aussi possible de faire suivre la determination des suites 
de droite d'une procedure de raccordement de droites. 

40 [0083] Ceci est represente sur la figure 15. Une solution peut consister a valider la droite D1 jusqu'a I'abscisse du 
premier point appartenant a la droite D2 et a recalculer la droite D2 a partir de la nouvelle valeur de grandeur dependante 
y ainsi determinee. En conservant le dernier point appartenant a D2, on determine ainsi une nouvelle droite D'2 repre- 
sentee en pointilie sur la figure 1 5. Une procedure analogue peut etre appliqu6e pour substituer la droite D'3 a la droite 
D3. On obtient ainsi un ensemble form6 de plusieurs droites de regression forme, dans I'exemple demerit, par les droites 

45 D1 , D'2, D'3. Cet ensemble constitue une approximation des grandeurs de mesure en reduisant une erreur maximale 
entre les grandeurs de mesure et Pensemble de droites. Une variante plus adapt6e a cette reduction consiste a faire 
que ce soient les droites correspondant aux plus fortes erreurs qui imposent leurs points limites comme points limites 
aux droites correspondant a des erreurs plus faibles. 

50 Dispositif cour une mise en oeuvre de la premiere phase 

[0084] Un schema d'un dispositif pour generer une fonction d'approximation selon invention est represente sur la 
figure 8. Le dispositif 5 recoit des couples (Xj, Yj) associant la grandeur dependante Yj a la grandeur independante X { . 
Les couples (X if Yj) entrent dans les premiers moyens 10 pour determiner et coder la fonction lineaire de regression 
55 formant une approximation des couples de mesure. Les codes specifiques ainsi determines sont transmis (connexion 
9) aux seconds moyens 17 qui determinent des seconds couples (X A , Y' A ) d'apres les codes specifiques a la fonction 
d'approximation. 

[0085] On distingue d'une part la seconde variante qui determine toutes les droites de regression passant par chaque 
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combinaison realisable de triplets avec la suite de premiers couples, et les autres variantes pour lesquelles a chaque 
droite de regression determinee (associee a son erreur E T ), on examine si les points additionnels restants fournissent 
bien des erreurs additionnelles interieures ou egales a I'erreur de triplet E T . 
[0086] La figure 9 represente un dispositif adapts a la seconde variante. II comprend : 

une memoire 12c MEM qui stocke notamment tous les points appartenant a I'ensemble de points a traiter. Ces 
points sont representes par leurs coordonnees (x, y), et eventuellement leurs poids W ou leurs inverses 1/W, 
une unite de calcul 13c COMPUT qui calcule, pour chaque triplet s6lectionn6, la droite de regression adapt6e a 
chaque triplet, c'est-a-dire les codes p, q de la droite et I'erreur E T associ6e au triplet. 

[0087] De plus, un controleur 11c CONTR permet d'effectuer la gestion des operations et d'adresser des nouveaux 
triplets en effectuant les lectures/ecritures de la memoire 1 2c et le chargement de I'unite de calcul 1 3c par de nouveaux 
triplets. La selection de la droite de regression qui est a conserver, c'est-a-dire pour cette variante ceile deiivrant la 
plus grande erreur de triplet, est effectuee par I'unite de calcul 13c. 

[0088] La figure 10 correspond au cas des autres variantes pour lesquelles pour chaque droite courante de regres- 
sion on examine si les autres points additionnels de la suite deiivrent une erreur inferieure a I'erreur de triplet. 
[0089] Pour cela, les premiers moyens 10 comprennent une unite de calcul 13c et une unite de comparaison 14c 
formant des moyens de calcul 1 9c. 

[0090] L'unite de calcul 1 3c, transmet les codes p, q de la droite de regression du triplet courant a I'unite de compa- 
raison 14c COMPAR qui determine si les points additionnels de I'ensemble de points generent une erreur plus faible 
avec cette droite de regression que celle generee par les points du triplet courant. Pour cela, I'unite de comparaison 
14c effectue le test : 

Si la droite de regression courante n'est pas accepted (test positif), un autre triplet est selectionn6 et une procedure 
analogue est a nouveau effectuee. Si le test est n6gatif pour tous les points additionnels, la droite de regression est 
acceptee et ses parametres sont charges par I'unite de calcul 13c dans la memoire 12c. 

[0091] Pour les variantes qui operent a partir d'une paire de points, I'unite de calcul 13c est programm6e pour de- 
terminer les points intermediates et pour former les triplets issus des paires de points, ceux-ci etant eux-m§mes even- 
tuellement issus d'enveloppes de points. Dans ce dernier cas, I'unite de calcul 13c est egalement programmee pour 
determiner les enveloppes. L'unite de comparaison 14c effectue ensuite les comparaisons des erreurs se rapportant 
aux differents triplets. 

[0092] Pour effectuer la gestion des operations et adresser des nouveaux triplets (lorsque des triplets examines 
n'ont pas donne satisfaction), le dispositif comprend, comme precedemment, un controleur 11c qui gere : 

les lectures/ecritures de la memoire 12c, 

le chargement de l'unite de calcul 1 3c par de nouveaux triplets, 

I'appel vers l'unite de comparaison 14c de tous les points additionnels a examiner successivement. 

[0093] II est possible de mettre en oeuvre la premiere phase de la methode a Paide de moyens de comparaison 14c 
(figure 12) ayant une structure neuronale. Consid6rons le cas le plus complet ou il existe des coefficients de pond6ration 
associes a chaque point P v Lorsque ces coefficients n'existent pas il suffit de leur donner une valeur unite dans les 
explications qui vont suivre. A titre d'exemple, on peut stocker dans la memoire 12c (figures 9, 10) les inverses des 
coefficients Wj affectes individuellement aux points P { . L'unite de calcul 13c determine selon les equations (2) et (3), 
pour un triplet donn6, les codes : -p, -q, E T et - E T . 

[0094] La condition a tester dans I'unite 14c pour un point additionnel P m a tester, de parametres X m , Y m , W m , est : 

Y m-PXm-(E T W m )-q>0 

ou 

Y m -P X m + <E T /W m )-q<0. 

[0095] Ces tests sont ais6ment realisables par une unite 14c ayant une organisation neuronale. En effet, ce test est 
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conslitue de fonctions lineaires, de fonctions de seuillage et de fonctions logiques aisement realisees dans un reseau 
de neurones. 

[0096] line telle unite neuronale de comparaison est representee sur la figure 12. Elle comprend trois neurones N1 , 
N2, N3. Les neurones N1 et N2 recoivent des donnees E1, E2, E3, E4. Le neurone N3 recoit les sorties des neurones 
N1 et N2. Chacune des entrees de ces neurones est aff ectee d'un coefficient synaptique Ci selon la technique connue 
mise en oeuvre dans les reseaux de neurones. Cette technique est par exemple decrite dans R.P. LIPPMANN, "An 
introduction to computing with neural nets" IEEE ASSP Magazine, April 1 987, pp. 4 to 22. 
[0097] Pour effectuer les tests precedents, on programme les neurones N1 et N2 selon le tableau I : 



TABLEAU I 



Neurone 


Ei 


c, 


E 2 


c 2 


E 3 


c 3 


E 4 


c 4 


1 




1 




-P 




-Et 


1 


-q 


2 


Y m 


1 


Xn, 


-P 


1/W m 




1 


•q 



[0098] Tous les coefficients synaptiques du neurone N3 sont egaux a 1 . 

[0099] Les neurones Nt et N2 calculent chacun un potentiel neuronal s tel que : s = I ( . C r E ( 

[0100] Ace potentiel s, il faut appliquer une fonction d'activation A. Pour les neurones N1 etN3, lafonctiond'activation 

At est : 

A 1 (s) = 1 sis >0 

A n (s) = 0 si s < 0. 
[0101] Pour le neurone N2, la fonction d'activation Ag est : 

A 2 (s) = 1 si s < 0 

A 2 (s) = 0 si s t 0. 

[0102] La sortie du neurone N3 vaut 1 si la condition a tester est satisfaite et vaut 0 dans le cas inverse. 
[0103] On observe sur le tableau I que les donnees -p, -q, -E T , + E T sont des codes qui apparaissent comme coef- 
ficients synaptiques dans le schema de la figure 12. Les donnees Y m , X m , 1/W m sont des donnees qui apparaissent 
sur les entrees du meme schema. Lorsque tous les points additionnels ont et£ testes et qu'une droite de regression 
a et6 s6lectionn6e les codes p, q, W charges dans la m6moire 12c sont alors utilises pour mettre en oeuvre la seconde 
phase de la methode (figure 11). 

[0104] L'interet presente par la realisation neuronale decrite est de pouvoir paralieiiser les diff6rentes operations a 
mettre en oeuvre selon les variantes deja decrites. Le fonctionnement d'une telle realisation neuronale est alors tres 
rapide. 

[0105] Lorsque les droites de regression ont et6 determinees, la seconde phase consistant a calculer des seconds 
couples (X A , Y' A ) de valeurs des grandeurs est alors mise en oeuvre dans des moyens de decodage 17 (figure 11). 
Les codes des droites sont charges dans une memoire 12a qui au cours de la seconde phase est adressee, par un 
controleur 11a, pour fournir les codes des droites de regression adressees. La memoire 12a, organisee en lignes par 
exemple, contient pour chaque droite de regression les parametres p, q et x L ou x L est la limite supeheure d'abscisse 
pour laquelle chaque droite de regression est d6finie. 
[0106] La memoire 12a contient ainsi un tableau de parametres 

x LvPc»qo 

X LnV Pm-V 

correspondant aux m droites de regression stockees. 
[0107] Les moyens de decodage 1 7 comprennent : 
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le contrdleur 11 a , 
la memoire 12 a , 

- et une unite de decodage 13 a . Celle-ci recoit unevaleurde requete X A de la grandeur independante pour laquelle 
on desire obtenir un tesultat Y' A issu de la fonction d'approximation. Pour cela le controleur 11 a presente succes- 
5 sivement certaines lignes de codes de la memoire 12 a et, pour chaque ligne de codes, I'unite de decodage 13 a 

teste le code x L pour determiner si la valeur X A de la grandeur d'entree peut etre traitee par !a droite de regression 
codee par ladite ligne. Si cela n'est pas le cas, I'unite 13 a requiert du controleur 11 a une autre ligne de codes. 
Lorsque cette condition est verifiee, I'unite 13 a effectue le calcul : 

10 Y' A = p n .X A+ q n 



si la ligne d'ordre n est la ligne adequate. 

is [0108] Comme cela vient d'etre decrit, les premiers moyens 10 peuvent constituer un dispositif de codage et les 
seconds moyens 17 peuvent constituer un dispositif de decodage, les deux dispositif s etant separes. lis peuvent en 
effet etre distants Tun de Pautre, les operations de codage etant separees des operations de decodage. Mais les deux 
moyens 1 0 et 1 7 peuvent aussi former un dispositif de codage/decodage unique. Dans ce cas certains organes peuvent 
etre affectes successivement a plusieurs operations. Notamment : 

20 

il peut s'agir d'une memoire 12 unique qui peut stocker les parametres des couples de valeurs et les codes des 

droites de regression determinees. Les unites 12a et 12c torment alors la memoire 12, 

il peut s'agir d'un controleur 11 unique qui regroupe les fonctions des controleurs 11c et 11a, 

I'unite de calcul 13c et I'unite de decodage 13a peuvent Sgalement former une seule unite 13. 

25 

[0109] II est egalement possible de realiser les moyens de decodage 17 (figure 11 ) selon une organisation a structure 
parallele. Une solution consiste alors a comparer, en parallele, la grandeur de requete X A a tous les codes x L . On 
obtient un ensemble de signaux T L repr6sentant les tesultats de ces tests de comparaison. On realise ensuite, en 
parallele, pour toutes les valeurs de L, un test qui combine deux a deux les signaux T L et T L+1 et qui indique si la 
30 grandeur de requete X A est dans I'intervalle d'indice L Ce test est verifie pour une unique valeur de L, notee n. On 
accede alors aux codes p n et q n correspondants, et on en d6duit la grandeur Y' A de la meme maniere que ptecedem- 
ment. 

[0110] Les premiers moyens 10 (de codage) associes aux seconds moyens 17 (de decodage) peuvent §tre utilises 
pour determiner la valeur d'une fonction d'approximation par au moins une droite de regression. Ces determinations 

35 peuvent etre effectuees pour des valeurs quelconques de la grandeur independante (a I'interieur de bomes pnSdefinies 
definissant I'etendue de Taction de chaque droite de regression). Cette methode evite de stocker inutilement des tables 
de valeurs, pour lesquelles toutes les valeurs ne seront pas exploitees. Selon la methode on ne determine que les 
seules valeurs necessaires a ('application. La methode selon invention presente I'interet de ne calculer que les valeurs 
necessaires. Le dispositif de I'invention peut etre utilise en combinaison avec un processeur neuronal pour calculer 

40 une approximation d'une fonction non lineaire par exemple une fonction sigmoide. Un processeur neuronal calcule 
des potentiels neuronaux issus de sommes de produits de coefficients synaptiques par des etats de neurones. La 
technique est developpee dans le document de R.P. LIPPMANN deja cite. Chaque potentiel neuronal doit etre soumis 
a Taction d'une fonction non lineaire. Pour cela, chaque potentiel neuronal constitue la grandeur independante decrite 
precedemment et le calcul de la fonction th(x) par exemple est effectu6 pour les seules valeurs de requetes utiles. La 

45 methode fournit des valeurs approximees, mais cette approximation peut etre obtenue avec une precision qui peut 
etre accrue selon le nombr-e de droites de regression calcutees et eventuellement les coefficients de ponderation W. 
La precision reste neanmoins tributaire de la precision des coordonnees des couples de depart. 
[0111] Une telle nrtethode est particulierement interessante pour les calculs de fonctions connues (telles que des 
fonctions mathematiques) ou de fonctions inconnues explicitement, par exemple une fonction representee par des 

so points de mesures, que I'on desire simplifier par une fonction lineaire de regression. 

[0112] Dans sa version neuronale, I'invention est interessante dans des applications neuronales car elle apporte 
non seulement un traitement homogene mais aussi une grande compacite de I'architecture materielle necessaire. 
[0113] II est possible de combiner les variantesde la methode, dans leurs mises en oeuvre mater ielles pour combiner 
leurs avantages. Ainsi par exemple, les combinaisons suivantes sontavantageuses : 

55 [0114] Dans le cas ou les poids existent, deux combinaisons apparaissent interessantes : 

Dans la premiere combinaison, on utilise tout d'abord la premiere variante basee sur des triplets de points. Si une 
solution n'est pas trouvee, on enchaine avec la deuxieme variante basee aussi sur des triplets de points. Ceci 
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permet d'obtenir les avantages suivants : 

Quelles que soient les donn§es, on obtient une solution. 

Dans le cas ou on exige qu'on obtienne toujours une solution exacte, cette combinaison permet d'obtenir la 
5 rapidite maximale. En effet, dans le cas general la solution est tournie directement par la premiere variante 

appelee, qui est celle lournissant le plus rapidement une solution exacte dans le cas ou les poids existent. La 
deuxieme variante est beaucoup plus lente mais elle permet de garantir que la combinaison fournit toujours 
une solution. 

io . Dans la deuxieme combinaison, on utilise tout d'abord la variante basee sur des paires de points. Si elle ne trouve 
pas une solution (qui peut etre approchee), on utilise la deuxieme variante basee sur des triplets de points. Ceci 
permet d'obtenir les avantages suivants : 

Quelles que soient les donnees, on obtient une solution (eventuellement approchee). 
15 . Cette combinaison est en moyenne plus rapide que la pr6c6dente. 

En contrepartie, elle ne fournit gSneVafement qu'une solution approchee. 

[0115] Dans le cas ou les poids n'existent pas, une combinaison interessante consiste a utiliser tout d'abord la va- 
riante basee sur I'enveloppe des points. Si elle ne trouve pas de solution, on utilise la deuxieme variante basee sur 
20 des triplets de points. Ceci permet d'obtenir les avantages suivants : 

Quelles que soient les donnees, on obtient une solution, 

Cette combinaison permet d'obtenir la meilleure rapidite moyenne. 

25 [0116] Lorsque la fonction lineaire de regression a ete determinee selon la methode decrite, il est possible de de- 
terminer non seulement des grandeurs Y' A ressortissant a ladite fonction mais egalement des grandeurs ressortissant 
a la derivee de ladite fonction. En effet, pour une valeur de requete X A de la grandeur independante, il suffit d'extraire 
des moyens de memorisation la valeur p ressortissant a la valeur X A pour connaTtre la valeur de cette d6riv6e corres- 
pondante. 

30 

Revendicatlons 

1. Dispositif (5) pour gene>er une fonction d'approximation fondee sur des premiers couples ((X-,, Y-,) a (X$, Y 6 )) de 
35 valeurs associant une grandeur dependante (Y., a Y 6 ) a une grandeur independante (X., a X$), et pour determiner 

des seconds couples (X A , Y' A ) de valeurs desdites grandeurs d'apres ladite fonction d'approximation caract6ris6 
en ce que le dispositif comprend : 

des premiers moyens (10) adaptes : 

40 

pour determiner iterativemant au moins une droite courante de regression D ayant pour equation D: y = 
p + q ou x et y sont des variables courantes, en rendant egales, en valeur absolue, des premieres erreurs, 
de signes altern6s, mesurees entre des premieres valeurs (Y 3 , Y 4 , Y 5 ) de la grandeur dependante pour 
trois couples (X 3 , Y 3 ) {X4, Y 4 ) (X 5 , Y 5 ) d'une suite desdits premiers couples, et respect ivement des se- 
45 condes valeurs (Y' 3> Y' 4 , Y' 5 ) de la grandeur dependante determin6es, d'apres ladite droite courante de 

regression, pour les memes valeurs (X 3 , X4, X5) de la grandeur independante, 

pour selectionner celle des droites courantes de regression qui deiivre I'approximation de tous les couples 
de ladite suite avec des erreurs minimales, 

et pour coder, a I'aide de codes specifiques p, q, la droite de regression seiectionnee, 

50 

et des seconds moyens adaptes (17) pour determiner lesdits seconds couples (X A> Y' A ) a I'aide desdits codes 
specifiques. 

2. Dispositif selon la revendication 1 caracterise en ce que lesdits premiers moyens (10) comprennent : 

55 

des moyens adaptes (12c) pour memoriser les grandeurs des premiers couples et lesdits codes specifiques 
pet q, 
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des moyens adaptes (1 3c) (1 9c) pour calculer successivement les codes specifiques p, q de droites courantes 
de regression et leurs premieres erreurs, egalisees en valeur absolue, respectives E T et pour selectionner au 
moins une droite courante de regression, 

et des moyens adaptes (11c) pour controler lesdits moyens (12c) pour memoriser et lesdits moyens (13c) 
pour calculer. 

Dispositif (5) selon les revendications 1 ou 2 caract6ris6 en ce que pour selectionner ladite droite de regression, 
lesdits premiers moyens comprennent des moyens (1 4c) adaptes : 

pour determiner pareillement des erreurs additionnelles pour les autres couples (X v (X2, Y 2 ) (Xq, Y 6 ) de 
valeurs de ladite suite, 

pour comparer les erreurs additionnelles avec lesdites premieres erreurs 6galis6es E T , 

et pour selectionner la droite courante de regression qui fournit des erreurs additionnelles inferieures ou egales 

en valeur absolue auxdites premieres erreurs egalisees E T . 

Dispositif selon la revendication 2 caracterise en ce que, pour selectionner la droite de regression, lesdits moyens 
adaptes pour calculer (13c) seiectionnent celle des droites courantes de regression qui d6livre des premieres 
erreurs egalisees E T maximales. 

Dispositif selon une des revendications 1 a 4 caracterise en ce que lesdits seconds moyens (17) comprennent : 

des moyens adaptes (12a) pour memoriser lesdits codes sp6cifiques p et q d'au moins une droite de regres- 
sion, 

des moyens adaptes (13a) pour selectionner, sur requete d'une valeur de grandeur independante X A , une 
droite de regression ressortissant a ladite valeur, et pour calculer, a partir des codes specifiques p et q recus 
des moyens pour memoriser, une valeur de grandeur dependante Y' A associee a la valeur de grandeur inde- 
pendante X A , 

et des moyens (11a) pour adresser aux moyens pour calculer, et sur leur demande, les codes sp6cifiques p 
et q des droites de regression. 

Dans un dispositif pour g6n6rer une fonction d'approximation fondee sur des premiers couples ((X,, Y t ) - (Xg, Y 6 )) 
de valeurs associant une grandeur dependante (Y-, - Y 6 ) a une grandeur independante (X 1 - Xg), m6thode pour 
generer ladite fonction d'approximation a partir des premiers couples f^, Y 1 )(Xe> Y 6 )de valeurs, et pour determiner 
des seconds couples (X A , Y' A ) de valeurs desdites grandeurs d'apres ladite fonction d'approximation, caracterisee 
en ce que, la methode comprend : 

une premiere phase : 

pour determiner iterativement au moins une droite courante de regression D ayant pour equation D: y = 
p + q ou x et y sont des variables courantes, en rendant egales, en valeur absolue, des premieres erreurs 
E PD , de signes alternes, mesurees entre des premieres valeurs (Y 3 , Y 4 , Y 5 ) de la grandeur dependante 
pour trois couples (X 3 , Y 3 ), (X^ Y 4 ) (X 5 , Y 5 ) d'une suite desdits premiers couples, et respectivement des 
secondes valeurs (Y' 3 , Y' 4 , Y' 5 ) de la grandeur dependante determinee, d'apres ladite droite courante de 
regression, pour les memes valeurs (X 3 , X4, X^ de la grandeur independante, 

pour selectionner celle des droites courantes de regression qui deiivre I'approximation de tous les couples 
de ladite suite avec des erreurs minimales, 

et pour coder la droite de regression seiectionn6e a I'aide de codes sp6crfiques p et q, 

et une seconde phase pour determiner lesdits seconds couples (X A , Y' A ) a I'aide desdits codes specifiques p 
et q. 

Methode selon la revendication 6 caracterisee en ce que la premiere phase comprend les etapes suivantes : 
A - selection (102) de trois couples (P3, P4, P5) de valeurs parmi ladite suite, 

B - calcul (104) de la droite courante de regression D et determination d'une erreur de triplet E^IEp^l s'y 
rapportant, 

C - selection (106) d'un couple addrtionnel (P1 , P2 ( P6), 
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D - calcul (108) d'une erreur additionnelle E Pm entre le couple additionnel et ladite fonction, 
E - et lorsque I E Pm I < £^110) pour le couple additionnel, la methode reprend a I'etape C avec un couple 
additionnel suivant, 

F - et lorsque I E Pm I > £ T pour au moins un couple additionnel (110), la methode reprend a i'etape A avec 
une selection d'un autre groupe de trois couples parmi ladite suite de couples, 

G - et lorsque I Ep^ I < E 7 pour tous les couples additionnels (112), la droite courante de regression est codee 
et stockee (114) en tant que fonction lineaire d'approximation. 

8. Methode selon la revendication 6 caracterisee en ce que la premiere phase comprend les Stapes suivantes : 
A - selection (102) de trois couples (P3, P4, P5) de valeurs parmi ladite suite, 

B - calcul (104) d'une droite courante de regression D et determination d'une erreur de triplet Ej=\E pt> \ sy 
rapportant, 

C - comparison (210) de I'erreur E T avec une erreur optimale E op initialised a une valeur strictement negative, 
D - et lorsque E r > E op , mise a jour (114) de I'erreur optimale E op en remplacant E op par E T et mise a jour des 

codes d'une droite optimale de regression D op en remplacant ceux-ci par les codes de la droite courante 

de regression D, 

E - puis retour (116) a I'etape A pour selectionner trois autres couples (102), 

F - et lorsque tous les triplets de couples de valeurs de la suite ont ete testes (116), les derniers codes de la 
droite optimale de regression D op constituent les codes de la fonction lineaire d'approximation (114). 

9. Methode selon la revendication 8 caracterisee en ce que la premiere phase comprend les etapes modifiees 
suivantes : 

A1 - modifie I'etape A en operant une selection (102a) de deux couples de valeurs appartenant a ladite suite, 
tel qu'il existe au moins un couple additionnel intermediate ayant une grandeur independante (X 1 - X$) 
comprise entre les grandeurs 

independantes dudit couple pour constituer au moins un triplet de couples de valeurs, 

A2 - modifie I'etape A premierement en determinant (102b) une droite annexe de regression qui contient les 
deux couples selectionnes et deuxiemement en determinant des secondes erreurs entre les grandeurs 
dependantes des couples intermediates possibles et ladite droite annexe de regression: 

et lorsque ces secondes erreurs sont toutes de meme signe, selection du couple intermediate four- 
nissant la plus grande seconde erreur, en valeur absolue, pour former un triplet de couples de valeurs 
fonm6 du couple intermediate et des deux couples selectionnes, 

et lorsque ces secondes erreurs sont de signes differents, reprise de la methode a I'etape A1 , 
B1 - I'etape B (104) est effectuee avec ledit triplet seiectionn6, 

C1 - modifie I'etape C (106) en seiectionnant (1 06a) un couple additionnel dont la grandeur independante n'est 
pas comprise entre la grandeur independante des deux couples selectionnes. 

10. Methode selon la revendication 9 modifiee en ce que, prealablement a I'etape A1 , lesdits couples de valeurs 6tant 
representes par des points (P1 - P6) dans un espace a deux dimensions, la premiere phase de la methode com- 
prend une etape (100a) de determination d'une enveloppe inferieure et/ou d'une enveloppe superieure reunissant 
les points les plus extremes de I'ensemble de points, la selection des deux couples de valeurs a I'etape A1 etant 
faite parmi les points adjacents appartenant a I'une ou ('autre enveloppe. 

11. Methode selon une des revendications 6 a 10 caracterisee en ce qu'a chaque valeur, prise dans la suite de valeurs 
de la grandeur independante, est associee un coefficient specifique de ponderation qui pondere lesdites premieres 
erreurs. 

12. Methode selon une des revendications 6 a 11 caracterisee en ce qu'elle genere une fonction d'approximation 
comprenant une suite de droites de regression raccordees entre elles. 



16 



EP 0 624 847 B1 



Patentanspruche 

1. Gerat (5) zum Erzeugen einer Annaherungsfunktion aut der Basis erster Wertpaare ((X 1( Y^ bis (Xg, Y 6 )), die eine 
abhangige Variable (Y 1 bis Y 6 ) mit einer unabhangigen Variablen [X^ bis Xq) verknupft, und zum Bestimmen zweiter 
Wertpaare (X A , Y' A ) der Variablen entsprechend der Annaherunqsf urktion, dadurch qekennzeichnet , daG das Gerat 
folgende Mittei enthalt: 

- erste Mittei (10): 

zum iterativen Bestimmen wenigstens einer aktuellen Regressionslinie D mit der Gleichung y = px + q, vyorin 
x und y aktuelle Variable sind, durch Angleichen erster Fehler mit wechselndem Vbrzeichen im Absolutwert, 
die zwischen ersten Werten (Y 3 , Y 4 , Y 5 ) der abhangigen Variablen fur drei Paare (X 3 , Y 3 )(X 4 , Y 4 )(X 5 , Y 5 ) einer 
Reihe der ersten Paare bzw. zweiter Werten (Y' 3 , Y 4 , Y' 5 ) der abhangigen Variablen gemessen sind, die 
entsprechend der aktuellen Regressionslinie fur dieselben Werte der unabhangigen Variablen bestimmt sind. 
zum Wahlen jener der aktuellen Regressionslinien, die die Annaherung aller Paare der Reihe mit Mindestfeh- 
lem erzeugt, und 

zum Codieren der gewahlten Regressionslinie mit Hilfe spezifischer Codes (p, q), und 

zweite Mittei (17) zum Bestimmen der zweiten Paare (X A , Y' A ) mit Hilfe der spezifischen Codes. 

2. Gerat nach Anspruch 1 , dadurch qekennzeichnet, daG das erste Mittei (10) folgende Elemente enthalt: 

Speichermittel (12c) zum Speichern der Werte der ersten Paare und der spezifischen Codes p und^qT* 
Rechenmittel (1 3c)(1 9c) zum aufeinanderfolgenden Berechnen der spezifischen Codes p, q aktueller Regres- 
sionslinien und ihrer betreffenden ersten Fehler mit Angleich im Absolutwert E T , und zum Wahlen zumindest 
einer aktuellen Regressionslinie, und 

Mittei (11c) zum Steuem der Speichermittel (12c) und der Rechenmittel (13c). 

3. Gerat (5) nach Anspruch 1 Oder 2, dadurch qekennzeichnet , daB zum Wahlen der Regressionslinie das erste Mittei 
Mittei (14c) enthalt: 

zum gleichartigen Bestimmen von Zusatzfehlern fur die anderen Wertpaare (X v Y^X^ YgKXe, Y 6 ) der Reihe, 
zum Vergleichen der Zusatzfehler mit den angeglichenen ersten Fehler E T , und 

zum Wahlen der aktuellen Regressionslinie, die Zusatzfehler ausgibt, die im Absolutwert kleiner als Oder gleich 
den angeglichenen ersten Fehlern E T sind. 

4. Gerat nach Anspruch 2, dadurch qekennzeichnet , daQ zum Wahlen der Regressionslinie die Rechenmittel (1 3c) 
jene der aktuellen Regressionslinien wahlen, die die maximalen angeglichenen ersten Fehler E T abgeben. 

5. Gerat nach einem oder mehreren der Anspruche 1 bis 4, dadurch qekennzeichnet , daB die zweiten Mittei folgende 
Elemente enthalten: 

Speichermittel (12a) zum Speichern der spezifischen Codes p und q zumindest einer Regressionslinie, 
Rechenmittel (13a) zum Wahlen einer zum Wert gehorenden Regressionslinie in Beantwortung eines Werts 
der unabhangigen Variablen X A , und zum Berechnen eines dem Wert der unabhangigen Variablen X A zuge- 
ordneten Werts der abhangigen Variablen Y' A , ausgehend von den spezifischen aus dem Speichermittel emp- 
fangenen Codes p und q, und 

Mittei (11a) zum Adressieren der spezifischen Codes p und q der Regressionslinie an die Rechenmittel bei 
Anforderung der spezifischen Codes p und q der Regressionslinien vom Mittei. 

6. Verfahren zum Erzeugen einer Annaherungsfunktion auf der Basis erster Wertpaare ((X 1f Y^ bis (X 6 , Y 6 )), die 
eine abhangige variable (Y1 bis Y 6 ) mit einer unabhangigen Nforiablen (X, bis X$) verknupft, und zum Bestimmen 
zweiter Wertpaare (X A , Y A ) der Variablen entsprechend der Annaherungsfunktion in einem Gerat zum Erzeugen 
einer Annaherungsfunktion auf der Basis erster Wertpaare ((X 1 , Y^ bis (X^, Y 6 )), die eine abhangige Variable (Y 1 
bis Y 6 ) mit einer unabhangigen Variablen (X 1 bis X^) verknupft, und zum Bestimmen zweiter Wertpaare (X A , Y A ) 
der Variablen entsprechend der Annaherungsfunktion, dadurch qekennzeichnet , daB das Verfahren folgende 
Schritte umfaBt: 

eine erste Phase: 
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zum iterativen Bestimmen wenigstens einer aktuellen Regressionslinie D mit der Gleichung D: y = px + 
q, worin x und y aktuelle Variable sind, mittels Angleichen erster Fehler mit wechselndem Vorzeichen im 
Absolutwert, wobei diese ersten Fehler zwischen ersten Werten (Y 3 , Y 4 , Y 5 ) der abhangigen Variablen 
fur drei Paare (X 3 , Y 3 ), (X4, Y 4 ), (X $I Y s ) einer Reihe der ersten Paare bzw zweiten Werten (Y' 3> Y' 4 , Y 5) 
der abhangigen Variablen gemessen sind, die entsprechend der aktuellen Regressionslinie fur dieselben 
Werte (X 3 , X 4 , X 5 ) der unabhangigen Variablen bestimmt sind, 

zum Wahlen jener der aktuellen Regressionslinien, die die Annaherung aller Paare der Reihe mit Min- 
destfehlern erzeugen, und 

zum Codieren der gewahlten Regressionslinie mit Hilfe spezifischer Codes p und q, und 
eine zweite Phase zum Bestimmen der zweiten Paare (X A , Y' A ) mit Hilfe der spezifischen Codes p und q. 

7. Verfahren nach Anspruch 6, dadurch gekennzeichnet , daB die erste Phase folgende Schritte umfaBt: 

A - Wahlen (102) von drei Wertpaaren (P3, P4, P5) aus der Reihe, 

B - Berechnen (104) der aktuellen Regressionslinie D und Bestimmen eines zugeordneten Tripelfehlers E T = 
lEppl. 

C - Wahlen (106) eines Zusatzpaares (P1, P2, P6), 

D - Berechnen (108) eines Zusatzfehlers E Pm zwischen dem Zusatzpaar und der Funktion, 

E - wenn IE PD I < E T (110) fur das Zusatzpaar (110) ist geht das Verfahren mit dem Schritt C eines folgenden 

Zusatzpaares weiter, und 

F - wenn IE PD I > E T fur wenigstens ein Zusatzpaar (110) ist, geht das Verfahren weiter mit dem Schritt A mit 
einer Auswahl aus einer anderen Gruppe von drei Paaren der Paarreihe, und G - wenn IE Pm I < E T fur alle 
Zusatzpaare (112) ist, wird die aktuelle Regressionslinie codiert und als Linearannaherungsfunktion gespei- 
chert(114). 

8. Verfahren nach Anspruch 6, dadurch gekennzeichnet, daB die erste Phase folgende Schritte umfaBt: 

A - Wahlen (102) von drei Wertpaaren (P3, P4, P5) aus der Reihe, 

B - Berechnen (104) der aktuellen Regressionslinie D und Bestimmen eines zugeordneten Tripelfehlers E T = 
IE PD I, 

C - Vergleichen (210) des Fehlers E T mit einem Optimalfehler E op mit einem strikt negativen Anfangswert, 
D - wenn E x > E op ist, Auffrischen (114) des Optimalfehlers E op durch den Austausch von E op gegen E T und 
Auffrischen der Codes einer optimalen Regressionslinie D op durch den Austausch dieser Codes gegen die 
Codes der Regressionslinie, 

E - anschlieBende Ruckkehr (116) nach Schritt A zum Wahlen von drei anderen Paaren, und 

F - wenn alle Tripel von Wertpaaren der Reihe uberpruft sind (116), bilden die letzten Codes der optimalen 

RegressionslinieD op die Codes der Linearannaherungsfunktion (114). 

9. Verfahren nach Anspruch 8, dadurch gekennzeichnet, daB die erste Phase folgende geanderte Schritte umfaBt: 

A1 - Andern des Schrittes A durch eine solche Wahl (102a) von zwei zur Reihe gehorenden Wertpaaren, daB 
sich zumindest ein zusatzliches Zwischenpaar mit einer unabhangigen variablen (X, - Xq) zwischen den un- 
abhangigen Variablen des Paares zur Bildung zumindest eines Tripelpaares befindet, 
A2 - Andern des Schritts A zunachst durch Bestimmen (102a) einer Zusatzregressionslinie, die das zweite 
gewahlte Paar enthat, und zweitens durch Bestimmen zweiter Fehler zwischen den abhangigen variablen der 
moglichen Zwischenpaare und der Zusatzregressionslinie: 

und, wenn diese zweiten Fehler alle dasselbe Vorzeichen haben, das Wahlen des Zwischenpaares zum 
Erzeugen des zweiten Fehlers mit dem groGten Absolutwert zur Bildung eines Tripelpaares der vom Zwi- 
schenpaar und den zwei gewahlten Paaren gebildeten Werte, 

und wenn diese zweiten Fehler verschiedene Vorzeichen haben, Fortsetzen des Vert ahrens im Schritt A1 , 
B1 - Schritt B (104) wird mit dem gewahlten Tripel durchgefuhrt, 

C1 - Andern des Schritts C (106) durch die Wahl eines Zusatzpaares, dessen unabhangige Variable nicht 
zwischen den unabhangigen Variablen der zwei gewahlten Paare liegt. 

10. Verfahren nach Anspruch 9, in dem Sinne geandert, daB vor dem Schritt A1 die Wertpaare'als Punkte (P1 - P6) 
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in einem zweidimensionalen Raum dargestellt werden, wobei die erste Phase des Verfahrens einen Schritt (100a) 
zum Bestimmen einer unteren Umhullenden und einer oberen Umhullenden umfa&t, die die AuGenpunkte der 
Punktgruppe verbinden, und die Wahl der zwei Wertpaare im Schritt A1 aus benachbarten zu einer der Umhul- 
lenden gehorenden Punkten getroffen wird. 

11. Veriahren nach einem oder mehreren der vorangehenden Anspruche 6 bis 10, dadurch gekennzeichnet, daft ein 
spezifischer Gewichtungskoelfizient, der das Gewicht der ersten Fehler bestimmt, jedem Wert in der Wertreihe 
der unabhangigen Variablen zugeordnet wird. 

1 2. Verfahren nach einem oder mehreren der vorangehenden Anspruche 6 bis 1 0 t dadurch gekennzeichnet, daB damit 
eine Annaherungsfunktion mit einer Reihe verknupfter Regressionslinien erzeugt wird. 



Claims 

1 . A device (5) for generating an approximation function based on first pairs ((Xj , Y-, ) to (Xe, Y 6 )) of values associating 
a dependent variable (Y, to Y 6 ) with an independent variable {X } to Xg), and for determining second pairs (X A , 
Y' A ) of values of said variables in accordance with said approximation function, characterised in that the device 
comprises: 

first means (10) adapted: 

to iteratrvely determine at least one current regression line D satisfying the equation D: y = px +q, where 
x and y are current variables, by making first errors of alternating sign equal in absolute value, which first 
errors have been measured between, respectively, first values (Y 3 , Y 4 , Y 5 ) of the dependent variable for 
three pairs (X 3 , Y 3 ) (X4, Y 4 ) (X 5 , Y 5 ) of a series of said first pairs, and second values (Y 3 , Y' 4 , Y' 5 ) of the 
dependent variable determined, in accordance with said current regression line, for the same values (X 3 , 
X4, X 5 ) of the independent variable, 

to select that one of the current regression lines which produces the approximation of all the pairs of said 
series with minimal errors, 

and to code the selected regression line with the aid of specific codes p, q, 
and second means (17) adapted to determine said second pairs (X A , Y' A ) with the aid of said specific codes. 

2. A device as claimed in Claim 1 , characterised in that said first means (10) comprise: 

storage means (12c) adapted to store the values of the first pairs and said specific codes, 
calculation means (1 3c) (1 9c) adapted to successively calculate the specific codes p, q of current regression 
lines and their respective first errors, equalised in absolute value (E T ), and to select at least one current re- 
gression line, 

and means (11c) adapted to control said storage means (1 2c) and said calculation means (1 3c). 

3. A device (5) as claimed in Claim 1 or 2, characterised in that for selecting said regression line said first means 
comprise means (1 4c) adapted: 

to likewise determine additional errors for the other pairs (X 1t Y,) (X 2 , Y 2 ) (Xg, Y 6 ) of values of said series, 
to compare the additional errors with said equalised first errors E T , 

and to select the current regression line which yields additional errors which, in absolute value, are smaller 
than or equal to said equalised first errors E T . 

4. A device as claimed in Claim 2, characterised in that for selecting the' regression line said calculation means (1 3c) 
are adapted to select that one of the current regression lines which yields maximal equalised first errors E T . 

5. A device as claimed in any one of the Claims 1 to 4, characterised in that said second means (17) comprise: 

means (12a) adapted to store said specific codes p and q of at least one regression line, 
means (13a) adapted to select, in response to a value of the independent variable X A , a regression line be- 
longing to said value, and to calculate, starting from the specific codes p and q received from the storage 
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means, a value of the dependent variable Y' A associated with the value of the independent variable X A , 
and means (11a) for addressing the specific codes p and q of the regression lines to the calculation means at 
their request. 

In a device for generating an approximation function based on first pairs ((X.,, Y-,) - (Xg, Y 6 )) of values associating 
a dependent variable (Y 1 to Y 6 ) with an independent variable (X 1 to Xg), a method of generating said approximation 
function based on first pairs (Xj, Y-,) (Xg, Y 6 ) of values, and for determining second pairs (X A , Y' A ) of values of 
said variables in accordance with said approximation function, characterised in that the method comprises: 

a first phase: 

for iteratively determining at least one current regression line D satisfying the equation D: y = px -tq, where 
x and y are current variables, by making first errors (E PD ) of alternating sign equal in absolute value, which 
first errors have been measured between, respectively, first values (Y 3 , Y 4 , Y 5 ) of the dependent variable 
for three pairs (X 3 , Y 3 ), (X4, Y 4 ), (X 5 , Y 5 ) of a series of said first pairs and second values (Y' 3 , Y' 4 , Y' 5 ) of 
the dependent variable determined, in accordance with said current regrssion line, for the same values 
(X 3 , X 4 , X 5 ) of the independent variable, 

for selecting that one of the current regression lines which produces the approximation of all the pairs of 
said series with minimal errors, 

and for coding the selected regression line with the aid of specific codes p and q, 
and a second phase for determining said second pairs (X A: Y' A ) with the aid of said specific codes p and q. 
A method as claimed in Claim 6, characterised in that the first phase includes the following steps: 
A - selecting (102) three pairs (P3, P4, P5) of values from said series, 

B - calculating (104) the current regression line D and determining an associated triplet error E T = IE PD , 
C - selecting (106) an additional pair (P1, P2, P6), 

D - calculating (108) an additional error Ep m between the additional pair and said function, 

E - if \EpJ < Fj{110) for the additional pair (110) the method proceeds with step C with a following additional 

pair, 

F - and if ) EPm \ > E T ior at least one additional pair (110) the method proceeds with step A with a selection of 
another group of three pairs from said series of pairs, 

G - and if IE Pm l £ E T ior all the additional pairs (112) the current regression line is coded and stored (114) as 
a linear approximation function. 

A method as claimed in Claim 6, characterised in that the first phase includes the following steps: 

A - selecting (102) three pairs (P3, P4, P5) of values from said series, 

B - calculating (104) the current regression line D and determining an associated triplet error E-^\E PD , 

C - comparing (210) the error E T with an optimum error E op having a strictly negative initial value, 

D - and if E T > E op updating (114) the optimum error E op by replacing E op by E T and updating the codes of an 

optimum regression line D op by replacing these codes by the codes of the current regression line D, 

E - subsequently returning (116) to step A to select three other pairs, 

F - and when all the triplets of pairs of values of the series have been checked (116) the last codes of the 
optimum regression line D op will form the codes of the linear approximation function (114). 

A method as claimed in Claim 8, characterised in that the first phase includes the following modified steps: 

A1 - modification of step A by effecting a selection (102a) of two pairs of values belonging to said series, such 
that there is at least one additional intermediate pair with an independent variable (X 1 - Xg) situated between 
the independent variables of said pair in order to form at least one triplet of pairs of values, 
A2 - modification of step A, firstly, by determining (102b) a subsidiary regression line which contains the two 
selected pairs and, secondly, by determining second errors between the dependent variables of the possible 
intermediate pairs and said subsidiary regression line: 

and, when these second errors all have the same sign, selecting the intermediate pair producing the 
second error having the largest absolute value in order to form a triplet of pairs of values formed by the 



20 



EP 0 624 847 B1 



intermediate pair and the two selected pairs, 

and, when these second errors have different signs, resuming the method in step A1 , 
B1 - step B (104) is effected with said selected triplet, 

C1 - modification of step C (106) by selecting an additional pair whose independent variable is not situated 
between the independent variables of the two selected pairs. 

10. A method as claimed in Claim 9, modified in that, prior to the step A1, said pairs of values being represented as 
points (P1 - P6) in a two-dimensional space, the first phase of the method includes a step (100a) to determine a 
lower envelope and/or an upper envelope which interconnect the outer points of the set of points, the selection of 
the two pairs of values in the step A1 being made from adjacent points belonging to one of the envelopes. 

11. A method as claimed in as claimed in any one of the Claims 6 to 10, characterised in that a specific weighting 
coefficient, which weights said first errors, is associated with each value in the series of values of the independent 
variable. 

12. A method as claimed in any one of the Claims 6 to 11 , characterised in that it generates an approximation function 
comprising a series of interlinked regression lines. 



21 



EP 0 624 847 B1 




22 



EP 0 624 847 B1 




FIG.3 
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FIG.9 




FIG.10 
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FIG.11 
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